Electron-Phonon Interactions in Correlated Systems: Adiabatic Expansion of 

Dynamical Mean Field Theory 
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We use the dynamical mean field theory to develop a systematic and computationally tractable 
method for studying electron-phonon interactions in systems with arbitrary electronic correlations. 
The method is formulated as an adiabatic expansion around the limit of static phonons. No spe- 
cific electronic groundstate is assumed. We derive an effective low- frequency phonon action whose 
coefficients are static local correlation functions of the underlying electron system. We identify the 
correct expansion parameters. At a critical electron-phonon interaction strength the system under- 
goes a transition to a polaronic state. We determine the location of this polaronic instability in the 
presence of electron-electron interactions, doping, and quantum lattice fluctuations and present the 
formalism needed for study of the electron self-energy and effective mass. 
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The development of a quantitative theory of electronic 
properties of "correlated materials" (i.e. those for which 
the local density approximation plus Boltzmann trans- 
port theory is inadequate) is an important goal of ma- 
terials science. An important feature of real materials is 
the electron-phonon interaction. The conventional the- 
ory of electron-phonon interactions in metals is due to 
Migdal and Eliashberg (ME)EI and is based on two as- 
sumptions: that the underlying electronic state is well 
described by Landau's Fermi liquid theory and that the 
typical phonon frequency luq is small compared to the 
electronic Fermi energy Ep , so that an expansion is pos- 
sible in the "adiabatic parameter" 7 = ujo/Ep. However, 
in many materials of current interest electron correla- 
tions are strong, so a Fermi liquid description may not 
be appropriate and the relevant expansion parameter is 
unclear. 

The introduction! and more recent improvements! of 
the dynamical mean field (DMF) method have opened 
an important avenue for progress, by showing how a 
good approximation to the correlation physics can be 
obtained from the solution of a quantum impurity prob- 
lem plus a self-consistency condition. Unfortunately the 
straightforward inclusion of the electron-phonon coupling 
in the DMF formalism presents a difficult technical prob- 
lem: the mismatch between the typical phonon frequency 
scale ujq < O.leV and electron energy scale t > leV 
renders conventional numerical approaches to the impu- 
rity problem prohibitively expensive, except in the "anti- 
adiabatic" limitQiJ 7 > 1 relevant to rather few materials. 

In this paper we show how to turn this apparent dif- 
ficulty to an advantage, by developing a systematic adi- 
abatic expansion of the DMF formalism about the limit 
7 = of static (classical) phonons which was shown to be 
easily tractable in Refs. [|[7| We derive an effective low- 
energy phonon action which, at leading order, reproduces 
ME theory, appropriately generalized to a ground state 
which is not necessarily a Fermi liquid. The simplifica- 



tions inherent in the DMF theory allow us to go beyond 
leading nontrivial order and to calculate the effect of dy- 
namic terms representing quantum lattice fluctuations. 
The vertices of the action are higher-order static local 
density correlation functions of the electron system in 
the absence of phonons. We identify the correctly renor- 
malized expansion parameters and show that the expan- 
sion breaks down in the vicinity of a polaronic instability, 
occurring at a critical value of the electron-phonon cou- 
pling. We determine the value of this critical coupling us- 
ing analytical and Quantum Monte Carlo methods. We 
also present the formalism needed to calculate the elec- 
tron self-energy. 

We write a general tight-binding based Hamiltonian as 
follows: H = H c \ + iJph + ifci-ph- The electronic part is 
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Here H oc represents electron-electron interactions, not 
explicitly written. The operator cj creates an electron 
with spin a on lattice site i. The chemical potential is fi 
and the mean density is n = (1/JV) ^2 icr c\ a Ci„, where N 
is the number of sites in the lattice. 

We model the phonons as quantum oscillators: 
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The operator Xi measures the ionic displacement at site 
i, M is the ion mass, and K is the spring constant. Wc 
have taken the oscillator frequency ujq — (K/M) 1 / 2 to be 
dispersionless (Einstein model) . To apply the method to 
more realistic situations one should interpret K and M 
as averages over the appropriate phonon bands or use,t 



extended DMF theory of Ref. 
can be trivially added and wil 
by the electron-phonon interaction. 

We take the electron-phonon interaction to be local: 



1. Anharmonic termsa'tB 
be seen to be generated 
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Here rii — ^ a c\ a Ci a> and (xi) — is the equilibrium 
phonon displacement for a uniform electron distribution. 

In DMF theory the properties of H may be obtained 
from the solution of an impurity model specified by the 
action S[c, c, x, a] = Sq[x] + S cc [c, c, a] + S\[c, c,x,a], with 



S o[ x ] = 7^J2 Xk ( K + MUJ D X -* 



(4) 



Si[c,c,x,a] — — y c nrT c na a n + g ^ c nr7 c n + kyCr x k . (5) 

Here S ee arises from H cc and is not explicitly written. 
The impurity electron and phonon are represented by 
c, c and x, respectively, and have been Fourier trans- 
formed according to x(t) — exp(—iuj k T)x k etc. Since 
the problem is local the fields depend on frequency 
but not on momentum. We use the compact notation 
x k = x(iu> k ), c n = c(iuj n ), with bosonic (fermionic) 
Matsubara frequencies to k = ZkirT (uj n = (2n + 1)itT) 
indexed by integers k (n). The local Green function 
Gioc is calculated from the impurity partition function 
Z[a] = J[dcdcdx] exp — S[c, c, x, a], 
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which defines the self-energy E. Notice that Z, G\ oc , and 
E are functionals of the mean field function a which is 
fixed by the condition that Gi oc [a] agrees with the mo- 
mentum integral of the full Green function, using the 
same self-energy, i.e. 



Gloc[ a ]n 



dek- 



fi - E[a]„ - e k 



(7) 



which depends on the lattice density of states p(e k ). 

The foregoing is general. To analyze the phonon prob- 
lem it is convenient to formally integrate out the electron 
fields 



exp(— S[x, a]) — / [dede] exp(— S[c, c, x, a]) 



(8) 



and work with an effective phonon action S[x, a] = 
+ Scc[ a ] + Si [a;, a]. We proceed by formally ex- 
panding S\ about the values x which extremize S. rThc 
crucial fact (of course already well-known to Migdala) is 
that the characteristic phonon frequency scale is small. 
Therefore, as we shall show, the frequency sums involv- 
ing phonon fields are dominated by frequencies of order 
luq = {K/M) 1 / 2 , thus Si may be evaluated by an expan- 
sion about the adiabatic limit. 

For small <?, x = 0; this is the conventional metal- 
lic state with no lattice distortion. As the coupling 
increases, eventually a state with x ^ becomes pre- 
ferred. This corresponds to a polaronic instability at 



which the ground state is fundamentally reconstructed 
and has been extensively discussedQ for models involving 
only clcctron-phonon interactions. For the rest of this 
paper we assume x — 0; i.e. we expand about the undis- 
torted ground state. We thus have: 



Si[x, a] = — tr In a 
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Y n [a\k 1 ,...,k n x kl ---Xk n , (9) 



fei, 



where the prime denotes the restriction to k\ + . . . + k n = 
0. The ^[0]*!,...,^ = T n [a\(iuj kl , . . . ,iuj kn ) are the con- 
nected n-point local density correlation functions of the 
electronic action. Their explicit form depends on the 
non-phonon physics contained in S cc . It is convenient to 
introduce a scale t on which electronic properties vary 
(this could be the bandwidth or some interaction scale) 
and to define the parameters 



7 



(K/M) 1 / 2 



t 



wo 
t ' 



A 



9 

Kt 



(10) 



We measure all frequencies and temperatures on the 
phonon frequency scale jt, writing uj n — Un/ijyt) and 
T = T/(jt), and we rescale the phonon fields to x k = 
x{iw k ) = (K/T) x l 2 xu so the free action becomes Sq[x] = 
\ J2k + &k)®-ki an d the interaction part becomes 



S± [x,a] = — tr In a — 



A n/2(«f 7 j«/2-l 



^ ' f „[a]fci,...,fc„£fci ■ ■ - Xk„ 
hi ,...,fe n 



(11) 



The parameter 7 now appears in the prefactor and 
in the frequency arguments of the rescaled vertex 
functions f „[a} kl ,...,k n = T„[a](iju kl , . . . ,ijQ kn ) = 
t n ~ 1 Tr n [a](iuj kl , . . . ,ioj kn ). Provided (as is usually the 
case) that the dominant contributions to T n come from 
frequencies of order t we may evaluate the r„ via a low- 
frequency expansion and may for most purposes set the 
external frequencies to zero and define coupling constants 
r n [o] := r„[a]o....,o- Notice that it is not necessary that 
the electron Green function has no structure on the scale 
of luq; merely that the dominant contribution to T n comes 
from the scale of t. With S given as in Eq. (|ll| ), E, G\ oc , 
and a may be computed from Eqs. (^), (|7j) as formal 
power series in 7. 

The complete action, up to terms of order 7 3 ^ 2 , is 



S[x, a] = -tr In a + - ^ x k D k 

k 

-lA 3 / 2 7 1 /2fi/2f 3[a] ikii 



'£~k<2*£' — k-L —k,2 



ki ,k2 



-A 2 7Tf 4 [a] ^2 Xk 1 Xk 2 Xk 3 x-k 1 -k 2 -k 3 - 

ki ,k 27 ks 



(12) 



2 



Eq. ( fL2[ ) can be regarded as an effective action on the low- 
energy scale 7*. The phonon Green function is D^ 1 = 1 + 
u)l-Xt 2 [a]k,-k ~ l-A/A c +a)2+A7a p |a) fe |+C'(7 2 ), where 
A c = r 2 [a]~ 1 is a critical interaction strength and a p is a 
damping parameter. The inclusion of static particle-hole 
diagrams r 2 in the phonon self energy reduces the phonon 
frequency scale but increases the electron-phonon inter- 
action strength. The renormalized expansion parameters 
are 

7 = 7(1 -A/AJ 1 / 2 , ^T^/V (13) 

Physical properties ~ A" l 7" (where m and n are the num- 
ber of phonon lines and loops, respectively, in the corre- 
sponding Feynman diagram) diverge like (1 — A/A c )"/ 2 ~ m 
at A — > A c (notice that m > n always), and the uni- 
form metallic-, state x = becomes unstable to local 
distortions.Enlil Here we consider only A < A c . 

The coefficients T n in fll^ ) are higher order static 
susceptibilities of the interacting electron problem and 
may be computed with the techniques presently avail- 
able. Here as an example we study the Holstein-Hubbard 
model with spin-^ electrons for which 

S CC = U dr[n T (T)-l/2]K(T)-l/2], (14) 
Jo 

where n a (T) = c a {T)c a {r). We specialize to a Bethe lat- 
tice for which p{e k ) = {At 2 - e 2 ) 1 ^^ 2 - e 2 )/{2nt 2 ) 
and the self-consistency equation (^) becomes a n = 
iu> n — t 2 G n . Here we focus on r 2 , but most considera- 
tions remain valid for T3 and T4,. Combining Eqs. (^) and 
( |TT| ) we can show that r 2 = tT((n — (n) cc ) 2 ) cc , the static 
local density-density correlation function, where the av- 
erage (• ■ -) cc is taken with respect to S cc . All calculations 
will be done in the paramagnetic (PM) state. 

In the weakly interacting limit u :=U / (2t) < lwe use 
standard perturbation theory. To second order in u we 
can replace the mean- field function a entering r 2 [a] by 
its noninteracting form ao = iu) n + /1 — t 2 a,Q ■ We find 

f 2 [a](u,/3,/i) = f 2 [a ](0, /?,//) - uf \ [a ] (0, 0, fi) 

+u 2 Qfi[a ](0,j9 > Ai)+CaC9,A*)) + 0(u 3 ), (15) 

where f 2[ao](0, (3, fx) = —(S/n) J^ 1 dxf(x — fj,)y/l — x 2 x 
is the noninteracting particle-hole bubble and G^{P, A*) is 
given as a sum of integrals which can be evaluated nu- 
merically. We set (3 = 2t/T and measure [i on the scale 
2t. Details will be given in a separate publication. 

In the strongly interacting limit u>lwe use Hartree- 
Fock (HF) theoryE 2 ] which approximates the local Green 
function by G n = |[(a„ - U/2)- 1 + (a„ + U/2)- 1 }. The 
resulting cubic equation for G n may be evaluated numer- 
ically and the result used to calculate r 2 . At half filling 
we find r 2 ~ j u ~ 3 for u — + 00. 



In order to study finite doping and the crossover from 
weak to strong electron-electron interactions we use the 
Quantum Monte Carlo (QMC) method as described in 
Ref. |^. The imaginary time interval is divided into 
L — (3/ At time slices. At each time slice an Ising spin 
is introduced to decouple the quartic U interaction. The 
partition function and all observables can then be sam- 
pled stochastically over all 2 L auxiliary spin configura- 
tions. Computer time grows with L 3 . The QMC algo- 
rithm has to be iterated self-consistently [using Eq. (0)] 
until a convergent solution for a is obtained. We found 
15-20 iterations to be sufficient. The starting guess for a 
was computed raising finite-temperature Iterated Pertur- 
bation TheoryB 

Fig. [l] shows the results of QMC simulations of T 2 as a 
function of u. Physical parameters are (3 = 6 and n = 0. 
The main graph shows the crossover from weak to strong 
coupling. In the insets we compare QMC (circles) to an- 
alytical results (lines), with excellent agreement. Upper 
inset: weak-coupling perturbation theory as discussed 
above, with f 2 [a ](0, 6, 0) = 0.745702 and C 2 (6,0) = 
—0.032789. Lower inset: logr 2 vs. log u which converges 
to the asymptotic solution logT 2 = —2 log 2 — 31ogu 
found above. 

In Fig. ^ we plot f 2 as a function of electron density n, 
for (3 = 6 and various u, computed by QMC with L = 32. 
For u = single-particle density fluctuations decrease 
when the system is doped away from half filling, as ex- 
pected from the Pauli principle. At u = 00, density fluc- 
tuations are forbidden at half filling, so T 2 (n = 1) = 0. 
The minimum at n = 1 develops as pairs (for n > 1) or 
holes (for n < 1) become the dominant density fluctua- 
tions at low doping. 

In general we expect A c also to depend on lattice vi- 
brations, i.e. non-zero 7. Integrating out cubic and quar- 
tic terms in Eq. (|l2| ) leads to an effective quadratic 
action S[x, a] = oX)fc^fc(l — •V^cOt) + ^t)^-k where 
K X ii) = r 2 + p7f 4 + C(7 2 ) at half filling and T = 0. 
The vertex f 4 [a ] = -16(1 - ^ 2 ) 3 / 2 (l - 6^ 2 )/(15vr) is 
negative at u = 0. A u — > 00 HF calculation gives 
F4 ~ — u~ 3 and no sign change. Thus T 2 (classical 
phonon self-energy) and T 4 (quantum lattice fluctua- 
tions) have competing effects: the former increases the 
coupling A — » A while the latter decreases A by renor- 
malizing A c — > A c (7). For superconductivity, a similar 
competition between "phonon dressing" and "vertex cor- 
rections" was observed in Ref. [l3|. 

We finally consider the electron self-energy £ = S co + 
S ph = Scj)/SG, where the Luttinger-Ward functional <j)[G\ 
is the sum of all vacuum-to-vacuum skeleton diagrams. 
The electron part S ee comes from <fi arising from See, 
the phonon part E ph from Si which may be expanded in 
powers of 7 as in Eq. (|l^) . To leading order we find 

^ = -^E^% k ^+ (7 2 )- (16) 
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In general S ph is of order 7 and its frequency dependence 
is on the scale of t so that it may be neglected compared 
to either the bare frequency dependence or S ee . However, 
if 6T/5G is singular at low frequencies, as in Fermi liq- 
uids, then dT,/du> may be of order unity. To investigate 
the possibility of such a singularity we note that 

^2[a]k,-k = — try^ G nrT G n+ k,a x 

(l + tT^2 Kin'k^n'a'Gn'+k^' ] , (17) 
n'a' / 

where the vertex A is given in terms of the particle-hole 
irreducible vertex A 1 via 



Following the usual Fermi liquid arguments we observe 
that A 1 is a smooth function of its arguments, so the re- 
quired singular behavior can only occur if we differentiate 
on one of the explicit G factors, leading to 
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with Aka = l+tTJ2 n , a , A^, G n , a ,G n , +k , a , and £P hre s a 
function which varies on the scale of t or U . If the ground 
state is a Fermi liquid then G na = — in sign(o> n )p(/x) + 
G lnc and we obtain for the phonon contribution to the 
mass enhancement 



m 
m 



l + \tp(fi)A 2 0a . 



(19) 
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We have calculated A and find it decreases rapidly as 
correlations increase, so the density-coupled electron- 
phonon interaction is "turned off" as the Mott insulator 
is approached. Details will be presented elsewhere. 

To summarize, we have combined a small-7 expan- 
sion with DMF theory to obtain a general formalism for 
studying electron-phonon effects in correlated electron 
systems. We have identified a local polaronic instabil- 
ity at A — > A c where physical properties diverge. We 
have shown that A c can be increased by electron-electron 
interactions and quantum fluctuations and decreased by 
doping. We have obtained formal expressions for the onc- 
phonon electron self-energy and effective mass. Future 
papers will apply the method to studies of conductivity, 
isotope effects, and dynamical consequences of electron- 
phonon interactions near a Mott transition. 

We thank S. Blawid and R. L. Greene for useful discus- 
sions. We acknowledge NSF Grant No. DMR00081075 
and the University of Maryland-Rutgers MRSEC for sup- 
port. 
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FIG. 1. Inverse critical coupling f 2 = A c 1 as a func- 
tion of electron-electron interaction u = U/(2t), at half 
filling and (3 = 6. Circles: L = 32 QMC results. 
Lines: weak-coupling perturbation theory (upper inset), 
strong-coupling Hartree-Fock theory (lower inset; notice 
log- log scale). 
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FIG. 2. QMC simulation of inverse critical coupling 
?2 = A^ 1 as a function of electron density n, for various 
values of electron-electron interaction u = U/(2t) and f3 = 6. 
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